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AN IMPROVED ENERGY ARGUMENT FOR THE 
HEGSELMANN-KRAUSE MODEL 

ANDERS MARTINSSON 


Abstract. We show that the freezing time of the d-dimensional Hegselmann- 
Krause model is O(n^) where n is the number of agents. This improves 
the best known upper bound whenever d > 2. 


1. Introduction 


The Hegselmann-Krause bounded confidence model, or simply the HK- 
model, is a simple model for opinion dynamics, first introduced in [T] and 
popularized in [2j- In this model, we consider n agents, indexed by integers in 
[n] = {1, 2,... , n}. Each agent i initially has the opinion xo(i), represented 
by a vector in V = for some d > 1. Two agents consider each others opin¬ 
ions reasonable if their Euclidean distance is at most a constant e, called the 
confidence radius. The agents update their opinions synchronously in dis¬ 
crete time steps by compromising with all opinions they consider reasonable. 
More precisely, for each t = 0,1,... we recursively define 
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where = {j € [n] : ||xt(i) — xt(j )||2 < e}. We will refer to (|1.1|) as the 

HK update rule. In this paper we will always assume that e = 1. 

Arguably, the most fundamental result about the HK-model is that, for 
any initial configuration, the system freezes after a finite number of time 
steps. That is, for sufficiently large T we have that xt = xt for any t > T. 
We will refer to the smallest such T as the freezing time of the system, and 
let Tfi{n) denote the maximal freezing time of any configuration of n agents 
with d-dimensional opinions. 

In [3] it is shown that, for any d, Td{n) = and is further conjectured 

that Td{n) grows polynomially in n. In dimension one, this was first shown 
in [3] who proved that Ti(n) = This has later been improved to 

Ti(n) = 0{n^) in [S], see also PE]. Polynomial freezing time in arbitrary 
dimension was shown in [S], which obtains the bound Tj_{n) = 0{n^^d?‘). As 
opinions will always be contained in the affine space spanned by the initial 
opinions, we may assume that d < n — 1. Hence, this implies a uniform 
upper bound of T(i{n) = 0(n^^) independent of d. In a recent paper [S], this 
was improved to Tj_{n) = 0{n^). 

The problem of finding lower bounds on Td{n) has received less attention. 
In [5| it was noted that it is possible to obtain freezing times of order 
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for any d > 2 by placing opinions equidistantly on a circle. More recently, 
[9] shows that a certain “dumbbell” configuration achieves freezing time of 
order n? also in d = 1. 

The aim of this paper is to prove the following upper bound on the freezing 
time. 

Theorem 1.1. The maximal freezing time for the n-agent HK model in any 
dimension is 0{n'^). 


2. Proof of Theorem 1.1 

For a given sequence {xt} € V”, we define the communication graph Gt 
as the graph with vertex set [n] and where i and j are connected by an edge 
if ||xt(z) — xt(j )||2 < 1. Note that all vertices in Gt have an edge going to 
themselves. We will here write i j to denote that i is adjacent to j in Gt- 
We further let Pt denote the transition matrix for a simple random walk on 
Gt- Hence, if we think of x^ as a n x d matrix, we can formulate the HK 
dynamics as 

(2.1) xt+i = Ptxj. 

For a configuration x = (x(l), x(2),... , x(n)) € V"" of n agents, we define 
its energy as 

n n 

( 2 . 2 ) T(x) = ^^min{||x(z) -x(j)|||,l} . 

i=l j=l 

Note that the energy of any configuration lies between 0 and Let {x^} be 
a sequence in V” which satisfies m- This energy function has the impor¬ 
tant property that £{^t) is non-increasing in t. The following inequality has 
played a central roll in obtaining the upper bounds on the high-dimensional 
freezing time in mi- 

Proposition 2.1. 

(2.3) £’(xt) - £:(xt+i) > 4||xt+i - xtll^. 

Proof. See Theorem 2 in [TU]. □ 

Here, we propose another way to estimate the energy decrement in a step 
in the HK-model. For a given state x and for any ordered pair (i,j) G [n]^, 
we say that {i,j) is active if ||x(f) — x(j )||2 < 1. We consequently define the 
active part of the energy of x as 

(2.4) Sactivefx) = ll^(*) “ ^(j)ll 2 

{ij) active 

Proposition 2.2. For each t > 0, let 

(2.5) \t = max{|A| : A 7 ^ 1 is an eigenvalue of Pt} . 

Then 

(2.6) £:(xt) - £:(x4+i) > (1 
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Proof. Let At denote the adjacency matrix of Gt, and let Dt denote its 
degree matrix, that is, the diagonal matrix whose (z,z):th element is given 
by the degree of i. Recall that every vertex in Gt has an edge to itself. 
Observe that Pt = At. We have 

^(^i) = Z] “ Mj)\\l + X] ^ 

= 2 Tr (xj {Dt - At) + X] 

i'/'tj 

where Tr(-) denotes trace. Here, we again interpret xt as an n x d matrix. 
Similarly, we have 

Sactivei^t) = “ ^t{j)\\l = 2 Tr (xj {Dt - At) Xt^ , 

ir^tj 

and 


^(xt+l)= Y1 ll^i+l(*) “^i+l0')ll2 + Y1 ^ 

= 2 Tr (xj^^ {Dt - At) Xi+i) + ^ 1 

i'/'tj 

= 2 Tf (x7 AtD^^ {Dt - At) D^^AtXt) + 1. 

i'/'tj 


Hence, it suffices to show that 

(2.7) Tr (x7 AtD^^ {Dt - At) D^^AtXt) < A? Tr (x^ {Dt - At) x*) . 

Let yt = D\^'^xt and Bt = D\^^PtD^ = D^ AtD^ . It is straight¬ 
forward to show that (j2.7l) simplifies to 

(2.8) Tr (y]Bt (/ - Bt) Btyt) < A? Tr (y7 {I - Bt) yt) • 

When d = 1, this inequality follows by standard linear algebra: write yt as a 
linear combination of eigenvectors of Bt and observe that Bt is a symmetric 
matrix which is similar to Pt. For the case when d > 2, let ei,..., denote 
the standard basis of V. We can rewrite (EID as 


(2.9) 


a a 

{yteY Bt {I - Bt) Btytet < Y [yJ{I 


2=1 


2=1 


Bt)ytei. 


By the one-dimensional case, we know that this inequality holds term-wise. 

□ 


Proposition 2.3. For any t >0, we have 

(2.10) |Ai| < 1 - 


diam(Gt) 
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where diam(Gt) denotes the graph diameter of Gt- If Gt is not connected, 
we interpret diam(Gi) as the largest diameter of any connected component 
ofGt. 

Proof. See for instance Corollary 13.24 in HD- Note that it suffices to con¬ 
sider the case where Gt is connected. □ 

Proof of Theorem \1.1\ . We call a time t = 0,1,... a merging time if two 
agents with different opinions at time t move to the same opinion at time 
t -|- 1. As merges are irreversible, there can at most be n — 1 such times. 

Assume that t < T is not a merging time. Then diam(Gt) > 2. Observe 
that for any i,j E [n], every second edge in a minimal path from i to j 
must have length at least hence SacUvei^t) = (diam(Gt)). Applying 
Proposition 12.21 it follows that the energy decrement in this step is It (;^), 
and can hence occur at most 0{n^) times. □ 

Remark 2.4. The energy argument presented here is optimal in the sense 
that there are configurations where the energy decrement is of order In 
particular, for the dumbbell in [9], this is the case during the first 0(n^) 
time steps until the communication graph changes. 
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